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CN ' We compare the numerical predictions for heavy quark production in high 

tZT* . energy hadron collisions of the conventional QCD parton model and the fc^- 

^^ i factorization approach (semihard theory). The total production cross sections, 

O I one-particle rapidity and pT distributions as well as two-particle correlations are 

c~| ' considered. The difference in the predictions of the two approaches is not very 

Mh, large, while the shapes of the distributions are slightly different. 
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1 Introduction 

The investigation of heavy quark production in high energy hadron colhsions is an 
important method for studying the quark-gluon structure of hadrons. Reahstic estimates 
of the cross section of heavy quark production as well as their correlations are necessary 
in order to plan experiments on existing and future accelerators as well as in cosmic ray 
physics. 

The description of hard interactions in hadron collisions within the framework of 
QCD is possible only with the help of some phenomenology, which reduces the hadron- 
hadron interaction to the parton-parton one via the formalism of the hadron structure 
functions. The cross sections of hard processes in hadron-hadron interactions can be 
written as the convolutions squared matrix elements of the sub-process calculated within 
the framework of QCD, with the parton distributions in the colliding hadrons. 

The most popular and technically simplest approach is the so-called QCD collinear 
approximation, or parton model (PM). In this model all particles involved are assumed to 
be on mass shell, carrying only longitudinal momenta, and the cross section is averaged 
over two transverse polarizations of the incident gluons. The virtualities g^ of the initial 
partons are taken into account only through their structure functions. The cross sections 
of QCD subprocess are calculated usually in the leading order (LO), as well as in the 
next to leading order (NLO) |l|, |^, |^, @, ^. The transverse momenta of the incident 
partons are neglected in the QCD matrix elements. This is the direct analogy of the 
Weizsaecker- Williams approximation in QED. It allows to describe quite reasonably the 
experimental data on the total cross sections and one-particle distributions of produced 
heavy flavours, however it can not reproduce, say, the azimuthal correlations of 
two heavy quarks, as well as the distributions of the total transverse momentum of 
heavy quarks pairs , which are determined by the transverse momenta of the incident 
partons. 

There is an attempt to incorporate the transverse momenta of the incident partons 
by a random shift of these momenta [kx kick) according to certain exponential 
distributions. This allows to describe quantitatively the two-particle correlations 0, 
but it creates the problems in the simultaneous description of one-particle longitudinal 
and transverse momentum distributions |^. 

Another possibility to incorporate the incident parton transverse momenta is referred 
to as /cy-factorization approach p, |10|, |Tl|, [12], |13[, or the theory of semihard interactions 



14| , p!5| , p!q , \[7\ , |18| [1^. Here the Feynman diagrams are calculated taking account of 
the virtualities and of all possible polarizations of the incident partons. In the small x 
domain there are no grounds to neglect the transverse momenta of the gluons, qiT and 
q2T, in comparison with the quark mass and transverse momenta, piT- Moreover, at very 
high energies and very high pi^ the main contribution to the cross sections comes from 
the region of Qit ~ Pit or g2T ~ Pit- The QCD matrix elements of the sub-processes 
are rather complicated in such an approach. We have calculated them in the LO. On 
the other hand, the multiple emission of soft gluons is included here. That is why the 
question arises as to which approach is more constructive. 



Most of the published papers on /cr-factorization have presented no numerical results 
or presented rather incomplete ones. Old sets of structure functions have been used, 
and, sometimes, the parton model results obtained with a particular set are compared 
with /cy-factorization results based on another set. 



In our previous paper [20| we have presented a comparison of results obtained with 



the help of fc-r-factorization and the parton model. The main goal of the paper pO 
was to demonstrate the differences in the qualitative and numerical predictions coming 
from the matrix elements. To simplify the calculations and to avoid various additional 
dependences we had used gluon distribution which had only a reasonable qualitative 
behaviour and a fixed value of as- 

The aim of this paper is to present a comparison between the results of the con- 
ventional parton model and the /cT^-factorization approach for the quantities which are 
measured experimentally. For this reason we use the realistic gluon distribution GRV94 
21|] compatible with the most recent data, see discussion in Ref. [Q. 

Below we shortly repeat the main formalism of the approaches used, discuss the val- 
ues of the parameters and present the numerical results on charm and beauty production 
obtained in the LO (and qualitatively in NLO) parton model and in the /cr-factorization 
approach. 



2 Conventional parton model approach 

The conventional parton model expression for the calculation of heavy quark hadropro- 



duction cross sections has the factorized form 23 



a{ah -^ QQ) =Y1 j dxidxjGa/i{xi, iJ,F)Gb/j{xj, iJ,F)d-{ij -^ QQ) , (1) 

where Ga/i{xi, fip) and Gb/j{xj, fip) are the structure functions of partons i and j in 
the colliding hadrons a and b, ^p is the factorization scale (i.e. virtualities of incident 
partons) and a{ij -^ QQ) is the cross section of the subprocess which is calculated in 
perturbative QCD. The latter cross section can be written as a sum of LO and NLO 
contributions, 

H^J^QQ) = ^|^(i;5'^(p)+4vra.(/i«)[i;«(p) + /;f^)(p)ln(/iVm|)]), (2) 

where Hr is the renormalization scale and /j ° as well as /j • and /j • depend only on 
the single variable 

_ 4m| ._ .^x 

(In the factor ln(/i^/mQ) we assume hr = fip following [l|]. In the case of different 
values of fj,R and fip, which is preferable for the description of the experimental data 
0, Eq. (2) becomes more complicated.) 



The expression (1) corresponds to the process shown schematically in Fig. la. The 
main contribution to the cross section at small x is known to come from gluons, i = j = 
9- 

Usually in the parton model the values 

Hf = IJ^R = rriQ (4) 

are used for the total cross sections and 



fiF = f^R = mT = -y/m^ + pI (5) 

for the one-particle distributions [^. However we calculate the total cross sections of 
heavy quark production as the integrals over their pt distrubutions, i.e. with scales (5). 

Both in the parton model and in the /cT^-factorization approach we take 

rric = 1.4 GeV, rrih = 4.6 GeV , (6) 



for the values of short- distance perturbative quark masses 

Another principal problem of the parton model is the collinear approximation. The 
transverse momenta of the incident partons, Qit and Qjt are assumed to be zero, and 
their virtualities are accounted for through the structure functions only; the cross sec- 
tion a{ij — > QQ) is assumed to be independent of these virtualities. Naturally, this 
approximation significantly simplifies the calculations. 

The conventional NLO parton model approach with collinear approximation works 
quite reasonably for one-particle diustributions and for the total cross sections; at the 
same time it is in serious disagreement with the data on heavy quark correlations (with- 
out kr kick introduction |7|). 



3 Heavy quark production in the /cT-factorization 
approach 

In the /cr-factorization approach the transverse momenta of the incident gluons in the 
small-x region result from the a^ In ^t diffusion in the gluon evolution. The diffusion 
is described by the function ip{x, q^) giving the gluon distribution at a fixed fraction of 
the longitudinal momentum of the initial hadron, x, and of the gluon virtuality, q^. At 
very low x that is to leading log(l/a;) accuracy, it is approximately determined |jl4| via 
the derivative of the usual structure function: 



^(.,,=) =4v^.»*^^, (7) 

Such a definition of (^^(x, q^) enables to treat correctly the effects arising from the gluon 
virtualities. 



Although generally (p is a function of three variables, x, qt and g^, the transverse 
momentum dependence is comparatively weak since g|> ~ —q^ for small x in LLA in 
agreement with g^-dependences of structure function. Note that due to QCD scaling 
violation the value ip{x, q^) for the realistic structure functions increases more fast with 
decreasing of x. Therefore at smaller x, larger qt becomes important in the numerical 
calculations. 

The exact expression of the qt gluon distribution can be obtained as a solution of 
the evolution equation which, contrary to the parton model case, is nonlinear due to 
interactions between the partons in the small x region. The calculations |^ of the qt 
gluon distribution in leading order using the BFKL theory [^ result in differences from 
our ip{x,q^) function given by Eq. (7) by only about 10-15%. 

Here we deal with the matrix element accounting for the gluon virtualities and 
polarizations. Since it is much more complicated than in the parton model we consider 
only the LO of the subprocess gg -^ QQ which gives the main contribution to the heavy 
quark production cross section at small x, see the diagrams in Fig. 2. The lower and 
upper ladder blocks present the two-dimensional gluon functions ^{xi, qf) and ip{x2, gf). 

Strictly speaking, Eq. (7) may be justified in the leading log(l/a;) limit only. To 
restore the unintegrated parton distribution /^(x, gy,//) (i.e. the probability to find a 
parton a with transverse momentum qt which initiates our hard process, with factoriza- 
tion scale /i) based on the conventional (integrated) parton density a{x, A^) we have to 
consider the DGLAP evolution^ 



da 



a. 



(91nA2 27r 



l-A 



Paa'{z)a'{-, X')dz - a{x, A^) ^ / Pa'a{z')dz' 



l-A 



(8) 



(here a(x, A^) denotes xg{x, A^) or xq{x, A^) and Paa'{z) are the splitting functions). 

The first term on the right-hand side of Eq. (8) describes the number of partons 6a 
emitted in the interval A^ < qf < A^ + 5A^, while the second (virtual) term refiects the 
fact that the parton a disappears after the splitting. 

The second contribution may be resummed to give the survival probability Ta that 
the parton a with transverse momentum qt remains untouched in the evolution up to 
the factorization scale 



Ta{qt,^) =exp 



^' asjpt) dpi ^ 

2 Z^ 

a' 



2vr pt 
Thus the unintegrated distribution fa{x,qt,fi) reads 



l-A 



Pa'a{z')dz' 



(9) 



fa{x,qt,fi) 



a 



^0(1-4-.) 



l-A 



X 



Paa'{z)a' ( -,gj ]dz 



Ta{qt,fi) 



(10) 



where the cut-off S = qt/ jj. is used both in Eqs. (9) and (10) |2^, ^ 



^For the g ^ gg splitting we need to insert a factor z' in the last integral of Eq. (8) to account for 
the identity of the produced gluons. 



In the leading log(l/x) (i.e. BFKL) limit the virtual DGLAP contribution is ne- 
glected. So Tfl = 1 and one comes back to Eq. (7) 

/r^-(x,g„;.) = ^^^ ,A^ = g,. (11) 

In the double log limit Eq. (10) can be written in the form 

/f ^^(x, qu /i) = jri—, fa(x, \')U\ /.)] ^_ , (12) 



ainA2 



A=(jt 



which was firstly proposed by ||30|. In this limit the derivative dTa/ d In X"^ cancels the 



second term of the r.h.s. of Eq. (8) (see [^ for a more detailed discussion). 

Finally, the probability fa{x,qt,^) is related to the BFKL function ip{x,qf) by 

^{x,q') = 4v^7rVa(x,gi,/i). (13) 

Note that due to a virtual DGLAP contribution the derivative da{x^\^)/d\^ can 
be negative for not small enough x values. This shortcoming of Eq. (11) is overcome 
partly in the case of Eq. (12). Unfortunately the cut-off A used in a conventional 
DGLAP computation does not depend on the scale /i. To obtain an integrated parton 
distributions it is enough to put any small A <^ 1^. 

On the other hand, in the survival probability Eq. (9) we have to use the true (within 
the leading log approximation) value A = gt//i. Thus for a rather large qt (of the order 
of /x) and x even the DDT form Eq. (12) is not precise enough. Only the expression (10) 
with the same cut-off A in a real DGLAP contribution and in survival probability (9) 
provides the positivity of the unintegrated probability fa{x,qt,f^) in the whole interval 
< X < 1. 

Of course, just by definition fa{x,qt,fi) = when the transverse momentum qt 
becomes larger than the factorization scale /i. 

In what follows we use the Sudakov decomposition for quark momenta pi^2 through 
the momenta of colliding hadrons p^ and pb {p'a = Pb — 0) and the transverse momenta 

Pl,2T- 

Pl,2 = Xi^2PB + yi,2PA + Pl,2T- (14) 

The differential cross section of heavy quark hadroproduction has the following form:^ 



dapp 1 



— J (f'qiTd^q2T5{qiT + q2T - Pit - P2t) 



dyldy^d^PiTd'^P2T (27r )8 (s)^ 

^There is a cancellation between the real and virtual soft gluon DL contributions in the DGLAP 
equation, written for the integrated partons (including all kt < A*) • The emission of a soft gluon with 
momentum fraction (1 — z) — > does not affect the x-distribution of parent partons. Thus the virtual 
and real contributions originated from 1/(1 — z) singularity of the splitting function P{z) cancel each 
other. 

On the contrary, in the unintegrated case the emission of soft gluon (with q'^ > kt) alters the trans- 
verse momentum of parent (i-channel) parton. Eq. (8) includes this effect through the derivative 
dTiklf^^)/dkl 

^ We put the argument of as to be equal to gluon virtuality, which is very close to the BLM 
scheme [pl|; (see also [|l7| ). 
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X ^^^{qlyMqlx)\MQQ\\ (15) 

Here s = 2paPb , <1i,2T are the gluon transverse momenta and y* 2 ^^^ ^^^ quark rapidities 
in the hadron-hadron c.m.s. frame, 

1/1 = ^ e^i, 2/2 = ^ eJ'S, y = yi + 2/2 (16) 

|MqqP is the square of the matrix element for the heavy quark pair hadroproduction. 
In LLA kinematic 

«?i ^ ypA + qiT, q2 ^ xpB + g2r- (17) 



so 



qf - -qfr, ql - -qlx- 



(The more accurate relations are qi = —j^-, qi = "iS but we are working in the 
kinematics where x, y ~ 0). 

The matrix element M is calculated in the Born approximation of QCD (see details 



m 



18, 20| without standard simplifications of the parton model. 



4 Total cross sections and one-particle distributions 

Eq. (15) enables us to calculate straightforwardly all distributions concerning one- 
particle or pair production. One-particle calculations as well as correlations between 
two produced heavy quarks can be easely obtained using, say, the VEGAS code p2| . 

However there exists a principal problem coming from the infrared region. Since 
the functions (p{x, g|) and ip{y, q\) are unknown at small values of g| and q\ we use the 
direct consequence of Eq. (7) p3| 



xG(x, q^) = xG{x, Ql) + ^-^ IJ^ ^(x, ql) dql (19) 



and rewrite the integrals in the Eq. (15) as 



/ (fqiTd^q2TS{qiT + q2T-piT-p2T)^^ '-Y-ip{ql,y)ip{ql,x)\MQQ\^ = 

J qi ?2 

= {AV2 7,'as{ml)f xG{x, Ql) yG{y, Ql) T\QI fi') L^ + (20) 




l-OO 

+ 4V2n^as{ml)xG{x,Ql)T{Ql,fi^) / dql^Siqir - Pit - P2t) 

Jq?, 



X 



Ql 




+ 4V2n^as{m'^T)yGiy,Ql)T{Ql,fi^) / dql^ 5{q2T - Pit - P2t) ^ 




X^^iqlx 

/■OO /"OO 

+ / d'^qiT / d'^q2T ^{qiT + q2T - Pit - P2t) x 
JQI JQl 

x^^^te^yVfe^x)|MQQ|^ 
Q'l Q2 

where the unintegrated gluon distributions are taken from Eqs. (9) and (10). 

The first contribution in Eq. (20), with averaging of the matrix element over the 
directions of the two-dimensional vectors qiT and q2Ti is exactly the same as the conven- 
tional LO parton model expression, with QCD scales /i^ = m^ and /i|. = Ql multiplied 
by the 'survival' probability T^{Ql, fi'^) not to have transverse momenta gi,2,T > Qo- 
The sum of the produced heavy quark momenta is exactly zero here. 

The next three terms contain the corrections to the parton model related to the gluon 
polarizations, virtualities and transverse momenta in the matrix element. The relative 
contribution of the corrections strongly depends on the initial energy. If it is not high 
enough, the first term in Eq. (20) dominates, and all results coincide practically (after 
accounting for Eq. (19)) with the conventional LO parton model predictions. In the case 
of very high energy the opposite situation takes place, the first term in Eq. (20) can be 
considered as a small correction and our results differ from the conventional ones. It is 
necessary to note that the absolute as well as the relative values of all terms in Eq. (20) 
strongly depend on the T-factor inclusion (i.e., when we use Eq. (10) or Eq. (7)). 

Before the numerical comparison it is necessary to note that the NLO parton model 
actually results only in a normalization factor in the case of one-particle distributions, 
the shapes of LO and LO-I-NLO distributions are almost the same, see 0, §, |], ^. This 
means that we can calculate the K-factor 

aiLO)+aiNLO) 
^= ^(LO) ' ^ ^ 



say, from the results for the total production cross sections, and restrict ourselves only 
to LO calculations oi pr, or rapidity distributions multiplying them by the K-factors. 

The numerical values of the K-factors depend ||35[ on the structure functions used. 



quark masses, QCD scales and initial energy, the dependence on the renormalization 
scale fiR being especially important. This seems to be evident, because the LO con- 
tribution is proportional to a^, whereas the NLO contribution is proportional to a;|. 
However the more important dependence at high energies, when small p values domi- 
nate, comes from the structure of Eq. (2). At p ^ the functions /j^-* and /^^^ have 
constant limits 0, /^^(p ^ 0) ^ 0.1 and /j])(p -^ 0) ^ -0.04, so due to Eq. (2) the K 
factor at high energies depends strongly on the ratio fi/rriQ. 

First of all let us present the role of the T-factors, Eq. (9). In Fig. 3 we show their 
values which were calculated as the ratios of the values of the first term of Eq. (20) to the 
same values calculated with Ta{qit, p) = 1 for the cases of charm and beauty production 
at -^/s = 14 TeV and p^ = s as functions of qu- The values of heavy quark transverse 
momenta were fixed at 20 GeV/c. In both cases the values of Ta{qit) are rather small 
at small qu and Ta{qu) -^ 1 at pt > qu- 

Now let us compare the numerical results predicted by the parton model and by the 
/cr-factorization approach. 

The energy dependences of the total cross sections of cc and bb pair production are 
presented in Fig. 4. As was mentioned, at comparatively small energies the first term in 
Eq. (20) dominates and the results of the fcr-factorization approach should be close to 
the LO parton model prediction. Actually the first results are even smaller due to the 
presence of the T-factor in Eq. (10). However the /cj^-factorization approach predicts a 
stronger energy dependence than the LO parton model both for cc and bb production. 
This can be explained by additional contributions appearing at very high energies in 
the fcy-factorization approach, see [pO[| 

The calculated values of one-particle pt distributions, da/dpT, in the /cT-factorization 
approach and in the LO parton model are presented in Fig. 5. In all cases the /c^- 
factorization predicts broader distributions. The average values of pt of the produced 
heavy quarks are rather different in these two approaches, as one can see from Table 1. 



Table 1 The average values of charm and beauty quark transverse momenta <Pt> 
(in GeV/c) in the fcr-factorization approach with /i2 = s and in the LO parton model. 





LO parton 


model 


A;T-factorization 


v^ 


cc 


bb 


cc 


bb 


14TeV 


1.78 


4.53 


2.23 


5.47 


1.8 TeV 


1.48 


3.96 


1.91 


4.54 



It seems to be very natural, because, contrary to the case of the LO parton model, 
a large p^ of one heavy quark can be compensated not only by the Pt of another heavy 
quark but also by the initial gluons (i.e. by hard gluon jet emission). 

The rapidity distributions of produced heavy quarks presented in Fig. 6 show that 
the main part of the difference between the fc-r-factorization approach and the LO parton 
model comes from the central region. 



5 Two-particle correlations 



We saw from the previous section that there is only a small difference in our results for 
the total cross sections and one-particle distributions obtained in the /cT-factorization 
and in the LO parton model. The predictions of the NLO parton model for these 
quantities differ from the LO parton model only by a normalization factor of 2-2.5 
IP, H, 1^, ^. So the difference between our predictions and the NLO parton model 
should be small. 

The calculations of two-particle correlations in different approaches are more infor- 
mative. The simplest quantity here is the distribution in the total transverse momentum 
of the produced heavy quark pair, Ppair- In LO, evidently, Ppair = Pit +P2T = Qit + 92T, 
and if qiT = q2T = 0, then da/dppair is a 5-function of zero. So the Ppair distributions give 
direct information about the transverse momentum distribution of the incident partons. 

It is clear that if q^T -C PiT, then the distributions in Ppair should be narrower in 
comparison with the one-particle pt distributions. In this case the Weizsaecker- Williams 
approximation should be valid and one can believe that the parton model reflects the 
real dynamics of the interaction. In the opposite case, q^r ~ PiT, the large transverse 
momentum of the produced heavy quark can be compensated not by the other quark, 
but by a high-p^ gluon. We have shown in our previous paper pO|, that about 70-80% 
of the total cross section of high-pT quark production at high energies originates from 
such processes, when the heavy quark propagator is close to the mass shell. 

We calculate the values of da/dppair for charm (a) and beauty (b) production in the 
fcr-factorization approach using the unintegrated gluon distribution Eqs. (9), (10) with 
scale values /i^ = s and /i^ = s/4 (only for ^/s = 14 Tev). Our results for pair production 
at different initial energies are shown by solid curves in Fig. 7. For the comparison we 
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present by dashed curves the one-particle pT-distributions taken from Fig. 5, obtained 
in the same fc^-factorization approach and with the same T-factor. As we put Qq = 1 
GeV^ in Eq. (20), we can not distinguish between the initial gluons with qx equal to, 
say, 0.1 GeV/c and 0.9 GeV/c, so our first bin in the da/dppair distribution has the 
width 2 GeV/c which explains some irregular behaviour of the solid curves at small pt- 
Naturally, all the solid and dashed curves are normalized equally at the same energy. 

At comparatively small energies, y/s = 39 GeV and even at y/s =630 GeV the 
distributions da/dppair are narrower than the one-particle distributions da/dpr- This 
means that the transverse momenta of the produced heavy quarks almost completely 
compensate each other. However the situation changes drastically with increasing of the 
initial energy. Starting from comparatively small pt, the difference between the curves 
decreases with energy. At ^/s = 14 TeV the distributions are similar both in the cases 
of cc and bb production. This means that the production mechanism changes in the 
discussed energy region. At ^i = 14 TeV the transverse momentum of the produced 
heavy quark is balanced more probably by one or several gluons (see also |2^), because 



the contribution with large virtuality in the quark propagator is more suppressed in 
comparison with the large virtuality in the gluon propagator. 

The discussed behaviour depends on the value of the scale /x^ in the T-factor, Eq. 
(9). The similar calculation at energy y/s = 14 TeV with //^ = s/4 is shown in Fig. 7 
for pair and single production by dash-dotted and dotted curves, respectively. Here the 
difference between these two curves is more significant and it becomes larger for lower 
energies. 

The distributions of the produced heavy quark pair as a function of the rapidity gap 
^y = Iuq ~ VqI between quarks are presented in Figs. 8. Here the difference between 
the LO PM and the fc^-- factorization predictions is not large again except for the region 
of very large Ay. 

Another interesting correlation is the distribution in the azimuthal angle 0, which is 
defined as the opening angle between the two produced heavy quarks, projected onto the 
plane perpendicular to the beam and defined as the xy-plane. In the LO parton model 
the sum of the produced heavy quark momenta projected onto this plane is exactly zero, 
and the angle between them is always 180". In the case of the NLO parton model the 
distribution is non-trivial [^ , however the predicted distribution (without including the 
k^ kick) is narrower in comparison with the fixed target data 0. 

The theoretical as well as experimental investigation of such distributions are very 
important to control our understanding of the considered processes. The problem is that 
in the case of one-particle inclusive distributions for heavy quark production in hadron 
collisions the sum of LO and NLO contributions of the parton model practically coincides 
P, |, ^ ^ with the LO contribution multiplied by K-factors. Therefore agreement with 
experimental data can be achieved for too small or too large NLO contribution by fitting 
one parameter, which can work as a normalization factor (say, the QCD scale). The 
deviation in azimuthal correlations from the trivial 6{(j) — tt) distribution comes from 
NLO correction to PM. However the standard NLO contributions are not large enought 
for the description of the data and only the comparatively large intrinsic transverse 
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momentum of incoming partons (Ajy-kick) allows to describe ^ the data. 

Preliminary results for the azimuthal correlations in the fc-r-factorization approach 
were considered in [|19|. The main difference in the information coming from da/dppair 



and da/d(f) distributions is due to the comparatively slow heavy quark. It gives a negli- 
gibly small contribution to the da/dppair in the first case, whereas in the second one each 
quark contributes to the distribution da/dcj) practically independently of its momentum, 
so all corrections coming from quark confinement, hadronization and resonance decay 
can be important. 

As was discussed above, the first contribution in Eq. (20) is the same as the con- 
ventional LO parton model in which the angle between the produced heavy quarks is 
always 180°. However the angle between two heavy hadrons can be slightly different 
from this value due to a hadronization processes. To take this into account we assume 
that in this first contribution the probability to find a hadron pair with azimuthal angle 
180° — is determined by the expression 

M<P) = ^== , (22) 

yPh + Pt 

where ph = 0.2 Gev/c is a transverse momentum in the azimuthal plane coming from 
the hadronization process. The other contributions of Eq. (20) result in a more or less 
broad distribution so we neglect their small modification due to hadronization. 

The /cT-factorization approach predictions for the azimuthal correlation of heavy 
quarks produced in pp collisions are presented in Fig. 9 and one can see that they 
change drastically when the initial energy increases from fixed target to the collider 
region. 



6 Conclusion 

We have compared the conventional LO Parton Model (PM) and the fcj^-factorization 
approach for heavy quark hadroproduction at collider energies using a realistic gluon 
(parton) distribution. Both the transverse momenta and rapidity distributions have 
been considered, as well as two-particle correlations, such as the distribution over the 
rapidity gap between two heavy quarks, their azimuthal correlations and distributions 
of the total transverse momentum of the produced heavy quark pair, Ppair- 

It has been shown in pO| that the contribution of the domain with strong qt ordering 



{Qi,2T "C rriT = Jm-g + Pt) coincides in fcr-factorization approach with the LO PM 
prediction. Besides this a very numerically large contribution appears at high energies in 
/cT-factorization approach in the region gi,2T > f^T- It kinematically relates to the events 
where the transverse momentum of heavy quark Q is balanced not by the momentum 
of antiquark Q but by the momentum of the nearest gluon. 

This configuration is associated with the NLO (or even NNLO, if both qi^2T > "^t) 
corrections in terms of the PM with fixed number of flavours, i.e. without the heavy 
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quarks in the evolution. Indeed, as was mentioned in [|I|, up to 80% of the whole NLO 
cross section originates from the events where the heavy quark transverse momentum 
is balanced by the nearest gluon jet. Thus the large "NLO" contribution, especially 
at large pt, is explained by the fact that the virtuality of the t-channel (or -u-channel) 
quark becomes small in the region qt — pr and the singularity of the quark propagator 
l/{p — q) —t^q) in the "hard" QCD matrix element, M{qiT,q2T,PiT,P2T), reveals itself. 

The double logarithmic Sudakov-type form factor T in the definition of unintegrated 
parton density (10) comprises an important part of the virtual loop NLO (with respect 
to the PM) corrections. Thus we demonstrate that /cT-factorization approach collects 
already at the LO the major part of the contributions which play the role of the NLO 
(and even NNLO) corrections to the conventional PM. Therefore we hope that a higher 
order (in as) correction to the fc^-- factorization would be rather small. 

Another advantage of this approach is that a non-zero transverse momentum of QQ- 
system {pTpair = Pit + P2T = Qit + 5'2t) is naturally achieved in the /cT-factorization. 
We have calculated the prpair distribution and compared it with the single quark pt 
spectrum. At low energies the typical values of pxpair are much lesser than heavy 
quark px in accordance with coUinear approximation. However for LHC energy both 
spectra become close to each other indicating that the transverse momentum of second 
heavy quark is relatively small. The typical value of this momentum {pTpair = A;r-kick) 
depends on the parton structure functions/densities. It increases with the initial energy 
(/c^-kick increases with the decreasing of the momentum fractions x, y carried by the 
incoming partons) and with the transverse momenta of heavy quarks, pt- Thus one gets 
a possibility to describe a non-trivial azimuthal correlation without introducing a large 
" phenomenological" intrinsic transverse momentum of the partons. 

It is necessary to note that the essential values of the parton transverse momenta 
qiT and g2T increase in our calculations with the growth of the value of p^**^ of detected 
6-quark. In the language of fc^-kick it means that the values of (fcf.) also increase. 

A more detail study of the heavy quark correlations in the /c^-factorization approach, 
the role of DL factor T(/c|,, ii^) and the value of scale n in T-factor will be published 
elsewhere. 
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Figure captions 

Fig. 1. Heavy quark production in liadron-liadron collision. The LO parton model 
corresponds to the case when qiT = q2T = 0. 

Fig. 2. Low order QCD diagrams for heavy quark production in pp (pp) collisions 
via gluon-gluon fusion (a-c). 

Fig. 3 The role of T-factor, Eq. (9) in the calculation of charm (solid curve) and 
beauty (dashed curve) production with pt = 20 GeV/c at ^/s = 14 TeV. 

Fig. 4. The total cross sections for charm and beauty hadroproduction in the 
/cr-factorization approach with unintegrated gluon distribution fg{x,qT,fi) given by 
Eq. (10), for /i^ values in Eq. (9) equal to s (solid curves), and s/4 (dash-dotted curves), 
and in the LO parton model (dashed curves). 

Fig. 5. pT-distributions of c-quarks (a) and 6-quarks (b) produced at different 
energies. Dashed curves are the results of the LO parton model. Solid curves are 
calculated with unintegrated gluon distribution /c,(a;, g-r, A*) given by Eq. (10), for /x^ 
values in Eq. (9) equal to s and dash-dotted curves are calculated for /x^ = s/4. 

Fig. 6. Rapidity distributions of c-quarks (a) and 6-quarks (b) produced at different 
energies. Dashed curves are the results of the LO parton model. Solid curves are 
calculated with the unintegrated gluon distribution fg{x,qT,fi) given by Eqs. (10), for 
/i^ values in Eq. (9) equal to s and dash-dotted curves are calculated for /x^ = s/4. 

Fig. 7. The distributions of the total transverse momentum Ppair for c-quarks (a) 
and 6-quarks (b) produced at different energies (solid curves), calculated with the un- 
integrated gluon distribution fg{x,qT,fJ') given by Eq. (10) and (10), for /x^ values in 
Eq. (9) equal to s. Dashed curves show the one-particle (single) pT-distributions with 
the same /i^ taken from Fig. 5. Dash-dotted and dotted curves are the same calculations 
for pair and single production at ^/s = 14 TeV with /i^ = s/4. 

Fig. 8. The calculated distributions of the rapidity difference between two c-quarks 
(a) or 6-quarks (b) produced at different energies in the kx factorisation approach, 
calculated with unintegrated gluon distribution fg{x,qT,fJ^) given by Eq. (10), for /i^ 
values in Eq. (9) equal to s (solid curves) and fj,"^ = s/4 (dash-dotted curves). Dashed 
curves show the LO PM predictioms. 

Fig. 9. The calculated azimuthal correlations of charm (a) and beauty (b) pair 
production for yU^ = s (solid curves) and fi"^ = s/4 (dash-dotted curves) at the energies 
equal to ^/s = 14 TeV, 1.8 TeV, 630 GeV and 39 GeV (the latter one only for charm 
production). 
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